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lattice QCD calculations. We examine the effects of mass dependent lattice artifacts by 
comparing two different formalisms for the heavy quarks. We take the continuum limit of 
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At a fixed lattice spacing, the mass of the charm quark is compared between quenched 
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limit of quenched QCD, we find m c (m c ) = 1.29(7)(13) GeV. No evidence was seen for 
unquenching. 
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1. Introduction 

The mass of the charm quark is a fundamental parameter of the Standard Model (SM), 
and yet its value is rather imprecise. The Particle Data Group Q quote 

1.15 < mlL(m charm ) < 1-35 GeV (1.1) 

This is to be contrasted with the more precise value of the mass of the D s : 

m Ds = 1.9695(5) GeV (1.2) 

from experiment. The problem is of course that quarks are confined hence their mass can 
never be directly measured. The imprecise value of the charm mass is just a reflection on 
how hard it is to solve QCD from first principles. 

The value of the charm mass is important for phenomenology. For instance the uncer- 
tainty in the charm quark mass is a big source of uncertainty in the production of charm 



- 1 - 



from DIS processes at HERA Q. Some models of physics beyond the standard model 
predict relations between various parameters such as the quark masses. The mass of the 
charm quark is also potentially important for understanding kaon decays |J. See the review 
article j|| for a comprehensive review of the mass of the charm quark. 

Lattice QCD can determine the hadron spectrum for a given quark mass. This can 
then be compared to experiment, and the quark mass tuned until the spectrum produced 
matches the experimental one. In practice the systematic uncertainties from the finite 
lattice spacing and too heavy sea quarks make this a non-trivial task in general. There 
have been many calculations of the mass of the charm quark from quenched QCD [||, |(| 
[7], ||, ^, 10 1, however there has not been an estimate of the error due ignoring the effect of 
virtual quark anti-quark pairs on the mass of the charm quark. 

What effect the sea quarks has on the mass of the charm quark is very difficult to 
estimate without simulating full QCD. However, using the running of the quark mass, in 
the quenched and Nf = 2 world, Mackenzie [11, 12] estimated that sea quarks could reduce 
the light quark masses by 10 — 20%. There have been claims that the light quark masses 



in unquenched QCD are significantly different to their values in quenched QCD. The CP- 
PACS collaboration [13] found that light quark masses with sea quarks were 25% lower 
than the quenched results. The recent computation, undertaken jointly by the HPQCD, 
MILC and UKQCD collaborations, of the strange quark mass is also significantly less 
than the result in quenched QCD [0, 15]. All the above calculations used perturbative 



renormalisation. 

Recently, there have been a number of two flavour unquenched lattice QCD calcula- 
tions, using Wilson or clover fermions, that have found that the strange quark mass is 
consistent with the value from quenched QCD [jO], 17, 18, FL9|. These new calculations 



use sea quarks with masses above a third of the strange quark mass, but do consistently 
use non-perturbative renormalisation techniques. The use of renormalisation factors, to 
two loop accuracy, by the HPQCD and UKQCD collaborations [S3] in the analysis of data 
from the lattice calculations that use improved staggered fermions moved the value of the 
strange quark mass closer to the quenched value, but still remained below it. The sea 
quarks used in the lattice calculations with Wilson or clover fermions [16, 17, 18, IS] are 
much heavier that those used by the HPQCD, MILC and UKQCD collaborations 14, 15]. 
The situation is not clear at the moment, but the introduction of sea quarks into a lattice 
QCD calculation could reduce the value of the strange quark mass by value between 0% 
and 10%. 

Only small differences have been found between the mass of the b quark in quenched 
and unquenched QCD p0| , 21, 22, 23]. The results for quark masses from the lattice have 
been reviewed by Lubicz 24], Wittig [25], and Rakow pq] . 

In this paper we make the first attempt to study the effect of sea quarks on the charm 
quark mass. Naively we would expect that the size of this effect would lie between that of 
the effect for strange and bottom quark masses. 

The mass of the charm quark is sizable in units of the lattice spacing that are compu- 
tationally feasible for unquenched calculations, hence the errors from the finite size of the 
lattice spacing are of great concern in this paper. There are a number of different lattice 
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K sea) 


Volume 


a 1 GeV ro = 0.5fm 


a 1 GeV ro = 0.55fm 


number of configurations 


(6.2,0) 


24 3 x 48 


2.913 


2.648 


216 


(6.0,0) 


16 3 x 48 


2.119 


1.926 


302 


(5.93,0) 


16 3 x 32 


1.860 


1.691 


278 


(5.2,0.1350) 


16 3 x 32 


1.876 


1.706 


395 



Table 1: Ensemble of gauge configurations. K soa = denotes a quenched ensemble. 



formalisms (recently reviewed by Kronfeld 27], Hashimoto and Onogi p3]), so there are a 



number of different ways of organising the calculations. 

Following the introduction, in section || we discuss the parameters of the lattice cal- 
culations. The different definitions of the quark masses are then outlined in section |3|, 
followed by section |], in which we describe the perturbative matching factors between the 
lattice data and the MS scheme. In section [| we discuss the methods used to interpolate 
from the quark masses at which the calculation was performed to the physical points. The 
final sections detail our results and conclusions. 

2. Details of the calculation 

The gauge fields were generated with Wilson's plaquette action and the quarks with the 
clover action, where the coefficient of the Sheikholeslami-Wohlert term, c$Wi has been 
determined non-perturbatively (NP). In this way the leading discretisation effects of the 
lattice are reduced from 0(a) to 0(a 2 ) for hadron masses, where a is the lattice spacing. 
Whilst this does not guarantee that lattice artifacts are smaller, the continuum limit is 
approached function of a . 

Hadron correlation functions were computed on three ensembles of gauge field config- 
urations in the quenched approximation, (3 = {6.2, 6.0, 5.93} and one ensemble of config- 
urations with the sea quarks, {(3 = 5.2, Ac sea = 0.135}. The values of the gauge coupling 
and quark mass of this ensemble were chosen so that the lattice spacing is matched to the 
coarsest quenched ensemble. The parameters of the lattice calculations are presented in 
table @. We will use the f3 value to distinguish each ensemble. The UKQCD collabora- 
tion has previously presented results and full details of the calculation for the light hadron 



spectrum and heavy-light spectrum and currents [3C] on the finest two of the three 
quenched ensembles, and the light hadron spectrum and currents on the matched quenched 
and sea quark ensembles We have already presented results for the heavy-light spec- 
trum on these four ensembles in [32]. Some of the charmonium mass spectrum from the 
unquenched data set has been reported in 33]. In this work we extend the analysis to the 
correlators necessary to define the mass of the charm quark. 

For the (3 = 6.0 and = 6.2 data sets, single or double exponential fits were made to 
the smeared source and local sink correlators. The gauge invariant smearing formalism of 



Boyle was used [34] with the parameters in [30|. For the (3 = 5.2 and (3 = 5.93 data sets, 



we fitted a variational multi-exponential fitting model to a 2 by 2 matrix of correlators 



made from a basis of local and fuzzed operators 3S] . 
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3. Definitions of the quark mass 



There are a number of different ways of calculating a bare lattice quark mass from the 
parameters in a lattice QCD calculation. The different definitions of the quark mass have 
different 0(a) effects. This is clearly seen in quenched calculations where a continuum limit 
is required for the two definitions to agree || Q . In this section we discuss various improved 
definitions of the quark mass that should have reduced lattice spacing dependence. The 
connection between the quark masses from the lattice and those in the continuum MS 
scheme is discussed in section ||. In the following we shall use lower cases for the quark 
masses and upper cases for the meson masses. 

One definition of the quark masses is from the Vector Ward identity 



1 1 1 

am ° = o 

Z \ K Kcrit 



(3.1) 



In the mass independent renormalisation scheme of the ALPHA collaboration 
the vector definition of the quark mass is 0(a) improved using 



3j, |9| 



arriQ 



am (l + b m am ) 



(3.2) 



The quark mass can also be defined from the Axial Ward Identity. This is often known 
as the PCAC mass. The axial current and pseudo-scalar density are defined as 



P(x) = 4li(x)j 5 1pj(x) 



(3.3) 



Although both the An and P operators depend on the flavor indices i and j, for simplicity 



we suppress the explicit dependence. The axial current can be improved according to |t0] 

A\{x) = A^x) + ac A d^P{x) (3.4) 
These currents are then renormalised as 



J 



R 



Zj(l + bj(amij))J l 



(3.5) 



where J is either A or P. The bare PCAC quark mass can then be defined in terms of 
correlation functions of the bare currents 



am. 



EMa\(x)p(o)} 

'" aC,l '~ 2£ x .(P I (x)P(0)) 



and the renormalised quark mass is given by 



am 



pcac,i] 



[1 + (b A - b p )am 0i ij] am pcaC)i j 



where the quark mass, m„ t ij, is given by 



am q ^ + am. 



2am 



(3.6) 



(3.7) 



(3.8) 
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for q either or pcac. 

The values of the improvement coefficients bj and cj are known to one loop in pertur- 
bation theory [^, [|l|. The improvement coefficients are also known non-perturbatively in 
quenched QCD for (3 > 6.0. The ca coefficient has recently been computed in unquenched 
QCD with 2 flavours of clover fermions We discuss the use of the non-perturbative 
improvement factors in section [|. 

To make the best use of the existing one loop results we use the tadpole improved 



formalism of Lepage and Mackenzie [43|. In this formalism the normalisation of the quark 
field changes, 

V2Kip -» V^uok^ (3.9) 

where uq is defined by 

u = (^tr[U^(x)}Y (3.10) 



Following Bhattacharya et al. jH, 45|, the expressions for the coefficients determined by 



Sint and Weisz can be re-written to form the tadpole improved expressions 

CA = -0.0952a s (3.11) 
bp = — (1 + 0.8763q s ) (3.12) 

Uq 

b A = — (1 + 0.8646a s ) (3.13) 
bm = — (-J-0.685a s ) (3.14) 

We also investigated other heavy quark formalisms that claim to have a smaller lattice 
spacing dependence. For heavy quarks with the improvement coefficients determined from 
one-loop perturbation theory the leading cut-off effects will be 0{a s am). In particular, for 
the dynamical ensemble, for which we cannot take the continuum limit, the lattice space is 
coarse (0(am) ~ 0.6). An effective field theory approach is the FNAL method |^] which 
supposes the dominance of 0((am) n ) cut-off effects over 0((ap) 2 ) effects. 

The lattice distorts the dispersion relation for a particle in the following way 

E 2 = Ml + ^ + (3.15) 
where M\ is the rest mass and M<i is the kinetic mass, defined by 

I d 2 E t . . 

w 2 = M lp=0 (3 - 16) 

An example of the energies of a typical data set fitted with both the continuum and FNAL 
dispersion relations is given in figure ||. 

In the Fermilab method it is the kinetic mass that is important for the dynamics of heavy- 
heavy and heavy-light bound states. The deviation of M\ from M 2 can give a measure of 
mass dependent cut-off effects. The quark mass, mi, is defined as 

ami = log(l + amo) (3-17) 
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aE 2 1-7 



1.6 



1.5 



FNAL dispersion relation 
continuum dispersion relation 




0.1 



0.2 



ap 



0.3 



Figure 1: Example of the energies for a typical data set (f3 = 5.2, kh = 0.1130 and kl = 0.1340), 
fitted using the continuum and FNAL dispersion relations. 



where too is defined in equation (|3.1[ ), Note the similarity of mi to the unrenormalised 
RGI mass in equation ( |3.2| ), to 0(a 2 ) at least, as 

(3.18) 



log(l + am) = am — - (am) 2 + - (am) 3 



and b m = — ^ at tree- level. This is the tree-level expression (in g 2 ) for the quark mass, to 
all orders in am. A perturbative definition of the kinetic quark mass is given by 



aTO2(aTOi) 



e™ 1 sinh(aTOi) 



(3.19) 



1 + sinh(aTOi) 

The quark masses, mi and TO2 can be used to get a perturbative definition of the hadron 
kinetic mass 

aM\ + (am-2 — ami) (3.20) 



aM 2 PT 



In this work we study equation 3.20, both at tree level and at one loop in perturbation 
theory. 

4. Perturbative matching 

We briefly describe the formalism required to extract the quark mass in the MS scheme at 



a specific reference scale. Most of the formalism is taken from |}47| , 36, 45]. We only use 
perturbative matching. Non-perturbative matching is discussed in section I 
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4.1 Quark mass renormalisation factors 

To extract the quark mass in MS we use the one loop matching factor 

CM = Z m {aix)Xm{(a) (4.1) 

where Z m {a\i) is the perturbative matching factor, X is tadpole improvement factor, and 
my i{a) is the bare lattice quark mass. 



1 _ a{p) s 

The QCDSF collaboration have published expressions for Cm as a function of the coefficient 



Z m (an) = 1 - -f^(81nOH - (C M - tad)) (4.2) 

4-7T 



of the clover term c sw [48]. 



C M = ^(12.952 + 7.738c SU) - 1.3804, ) (4.3) 

The numerical value of Cm = 25.758 for Csw = 1- At one loop it is consistent to use the 
one loop value for csw- Recently, the QCDSF collaboration have claimed substantial differ- 
ences between the renormalisation of the vector definition of the quark mass depending on 
whether the singlet or non-singlet estimate of the Z m factor is used []l6| . In this calculation 
there are no charm quarks in the sea, hence only the standard non-singlet renormalisa- 
tion factor is used for the vector quark mass. The ALPHA collaboration compute the 
connection between the quark mass on the lattice and the renormalisation group invariant 



mass [39| using a non-perturbative procedure. 

The tadpole improved |4^, value for Z m , based on X = 8n cr i t , uses tad = 10.66. 
For the non-tadpole improved case: X = 1 and tad = 0. 



At one loop, in some sense the scale \i in equation 4.1 is a free parameter. In principle 
no physical prediction should depend on the value of fi. The dependence on /i is reduced 
as the number of loops is increased. Reasonable choices for // lie in the range from \ ja to 
it J a. The "best guess scale" for the \i (called q*) , that attempts to minimise higher order 
corrections, can in principle be computed using the formalism described by Lepage and 



Mackenzie [43|. There has been a recent calculation of the q* for many of the perturbative 



expressions required for the PC AC mass [49] 



The connection between m MS and the PCAC mass is in |4~4 



in. 



MS 



The tadpole improved matching factors for the axial and pseudo-scalar operators are [44]: 
Z p (jj) = n (l + a s (J^ log(au) 2 - 1.328^ ) (4.5) 

and 

Z A = u (l -0.416a s ) (4.6) 
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The Fermilab group [ 50 1 have computed the connection between the lattice quark 
mass and the pole quark mass to all orders in m at one loop order in the expansion of the 
coupling. 

mi = mf ] + g 2 Z^l tanh mf ] (4.7) 

where 



Z& = ^\e- mf) coshmi 0) (^ (mi°VW(smh(mW)) -^(mi 0) )i?^(sinh(mW))) (4.8) 



,(o), 



,(o) 



(0)> 



,(0)^ 



The Zm\ factor is a function of the mass and the clover coefficient that can be reconstructed 
from the coefficients of the Chebyshev polynomials [50]. The functions Apy, Bpy, pa and 
pc are functions explicitly quoted in the Fermilab paper [50|. In the limit m\ — > 0, 



Z« = 0.143 + 0.0653c sw - 0.0116c 2 



mi 



^ g (m?f 



(4.9) 



This is the lattice part of the matching factor above. The Z^l factor contains the b m factor 
of Sint and Weisz [f£l}| . In the static (infinite mass) limit m^\oo) = 0.168<7 2 [51|. 
The one loop expression for the mi mass is 



m 



(i) 



m 2 {m [ i + g I m { {){l+g z Zm 



(4.10) 



where the function mi is defined in equation 3.19. The function Z M i is a function of m\ 

in the paper by the Fermilab group [50]. 

The tadpole improved definition of the m\ mass is 



m 



(i) 



+ 



M 



1 + M 



(4.11) 



1/12. 



where u\ 

The mi and mi mass definitions advocated by the Fermilab group are the lattice part 
of the matching. To convert the results to the MS scheme, the lattice results have the log 



term subtracted from equation 4J} and the results are multiplied by ■^fnal->ms(/ x ) - 



V FNAL— >MS 



00 



1 



3vr 



(4 + 61og()uo)) 



(4.12) 



The importance of subtracting the log term in equations of the form 4J; has been stressed 
by Groote and Shigemitsu [B2[. This is equivalent to the continuum matching factor used 



by Davies and Thacker [53 and so the matching factor agrees with that in equation 4.2 
In our early presentations |)4|, |55| we did not subtract the log term from equation 4.9 and 
only used the matching factor for the pole mass to the MS scheme. 

4.2 Discussion of non-perturbative matching 

So far we have only discussed the use of matching factors to one loop in perturbation 
theory. There are a number of elegant numerical formalisms that can compute renormal- 
isation factors non-perturbatively (see the reviews 

H, H II)- These methods promise 
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a non-perturbative matching with an accuracy limited by continuum calculations that are 
usually known to at least two loops. There is also a program of research into comput- 
ing perturbative matching factors, between the lattice and continuum QCD, to two loop 
accuracy [|5£fl. 

A particularly nice example of the power of non-perturbative matching was the com- 
putation of the mass of the charm quark in quenched QCD by Rolf and Sint 0. The 
consistent use of the non-perturbative factors from the ALPHA collaboration coupled with 
a controlled continuum extrapolation produced a very precise value for the charm mass in 
quenched QCD. 

There are now many non-perturbative estimates for matching factors from unquenched 
QCD with clover fermions |]60| , 61, 62, [l6|, |4^, |jf]. However, at the lattice spacing of 



our unquenched data it is not clear that non-perturbative renormalisation factors should 
automatically be used, unless the data is part of a consistent continuum extrapolation. 

In their calculation of the charm mass in quenched QCD, Rolf and Sint Q used im- 
proved coefficients and renormalisation factors determined non-perturbatively by the AL- 
PHA collaboration in quenched QCD. At non-zero lattice spacing there were significant 
differences between the vector and PCAC quark masses, that extrapolated to zero as the 
continuum limit was taken. Rolf and Sint performed their lattice calculations with lattice 
spacings in the range 0.1 to 0.05 fm. 

The unquenched data used in this work is at a fixed lattice spacing of 0.1 fm. At 
the moment it is computationally prohibitive, to do unquenched calculations with light 
sea quark masses and clover fermions, with a lattice spacing of 0.05 fm, using existing 
algorithms and computers [64, 65]. 

The non-perturbative estimates of improvement and matching factors can make the 
0(a 2 ) corrections to the quark sizable. There is an 0(a) ambiguity to the renormalisation 
condition used for the non-perturbative estimate of improvement coefficient or matching 
factor. Different conditions can produce different results at non-zero lattice spacing, but 
they will agree in the continuum limit. 

It has been observed that there is a large 0(a) ambiguity in the coefficient ca [44, 
45]. This induces an 0(a 2 ) ambiguity into the currents. This was noted because of the 



disagreement between the perturbative and non-perturbative estimate of ca at a lattice 
spacing of 0.1 fm in quenched QCD. 

The ALPHA collaboration have also found an 0(a) lattice spacing ambiguity in the 
coefficients bA — bp |6^| . Provided that a consistent definition of the coefficients is used 
at each lattice spacing, the continuum will be approached smoothly, and the different 
definitions should have the same continuum limit. A discussion of the effect of different 
determinations of the improvement coefficients can be found in pOj] . In particular, a naive 
comparison of the continuum limit of the quenched decay constant fx , with that obtained 
using the ALPHA determinations [67], using different determinations of the improvement 
coefficients suggests that this is indeed the case. 

In quenched QCD, the tadpole improved formalism has been extensively compared 
against the non-perturbative results [49]. A difference of 4% between the non-perturbative 
estimate and improved perturbative estimate was claimed. For our unquenched data at (3 = 
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5.2, we saw a 10% difference between the tadpole improved estimate of and the recent 
non-perturbative estimate [fj3]. Other groups have claimed to see systematic differences 
between using perturbative and non-perturbative renormalisation factors [16, 17, 19]. 

All current estimates of non-perturbative matching factors are done at leading order 
in the quark mass. At the lattice spacings we work at, the 0{am) terms are not small. 
In figure [2] the various renormalised quark masses are plotted. The coupling for f3 = 6.0 
is used. The renormalisation factors are expanded in both quark mass and coupling. The 
ALPHA formalism only treats the quark mass renormalisation to leading order. A large 
deviation of the FNAL renormalised quark masses from the masses renormalised using 
the ALPHA method would show that a one loop renormalisation factor with all orders in 
the quark mass could be closer to the continuum result, than using the ALPHA analysis, 
even with non-perturbative matching, at fixed lattice spacing. Of course the ALPHA 
renormalisation method is better than the FNAL renormalisation as the lattice spacing 
is reduced, but this may not be computationally feasible. There have been attempts to 
develop a non-perturbative definition of the FNAL formalism [68, 39 1. 



1.5 



am i 



0.5 



O FNAL ml 

□ FNAL m2 

O ALPHA numerical 

* ALPHA perturbative 




am, 



Figure 2: Renormalised group invariant quark mass as a function of the vector quark mass at /3 = 
6.0. The bursts are the FNAL renormalisation, the diamonds are the renormalised mass from the 
numerical calculations of ALPHA || , and the squares are the perturbatively renormalised masses 
from the perturbative expressions in the ALPHA formalism. 

A naive application of the ALPHA formulation at a fixed lattice spacing can be prob- 
lematic for heavy quarks. For example in figure [2| the renormalisation factors from ALPHA 
bend over at masses am ~ 0.6. These quark masses are too heavy for the use of 0(a) im- 
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provement at this lattice spacing. If calculations at quark masses larger than the charm 
mass are required, then they should be done at a finer lattice spacing. An alternative 



strategy is do a heavy quark interpolation in the continuum limit [7C]. 

Although Rolf and Sint have demonstrated that the ALPHA formulation can be used 
to compute the mass of the quark charm using data with lattice spacings at and finer than 
0.1 fm Q, the lattice spacing errors at 0.1 fm are not small. Given the high computational 
cost of reducing the lattice spacing errors in unquenched lattice QCD calculations, we 
feel that it is more useful to use tadpole improved perturbation theory to one loop and 
investigate the FNAL formalism for this data set. This should give a result closer to the 
continuum limit. As the lattice spacing is reduced in unquenched calculations, then the 
non-perturbative renormalisation will be crucial in producing results with high accuracy. 

4.3 Evolving the quark mass to the charm quark mass 

The matching of the quark mass in the lattice scheme to the quark mass in the MS scheme 
produces the mass at the scale jx, where fi is chosen, or guessed, to minimise the higher 
loop corrections to the perturbative matching factor. It is conventional [Q to evolve the 
quark masses from the scale /itoa standard reference scale of m c GeV. This is sometimes 
known as the scale invariant mass. 

The evolution is done using the solution of the renormalisation group equation 

^m ( "'V) = lm f m^ n f\fi) (4.13) 

given by 

^0) = c(a s (n)/ir) ^ 
m(/jLo) c(a s (//o)A0 

The anomalous dimension 7^/ is known to four loop order. We use the RunDec [71] 
mathematica package to do the evolution. 

4.4 Coupling prescription 

For the perturbative matching a choice of coupling is required, or equivalently a choice 
of Aqcd- We use the Ams computed on the same data set W^. This allows us to use a 
consistent coupling in all stages of the calculation. This is sometimes known as "horizontal 
matching" fTql - (The j3 = 5.95 result is used for the matched quenched data). These 
values are also partially quenched. We consistently use rif = 2(0) in all the perturbative 
expressions for the dynamical (quenched) data set. The couplings used are presented in 
table |2[ The coupling a j^g( ( z) was calculated at any scale using the standard four loop 



evolution equation [71]. 



5. Interpolations in quark mass 



We determined the heavy-light hadron spectrum [32] and the four definitions of the quark 
mass for each combination of heavy and light k value. We used the mass of the D s (1.9685 
GeV) meson to determine the mass of the charm quark. To avoid a large extrapolation 
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p 




a' 1 GeV 


A^MeV 


a s {a~ x ) 


a s (7ra -1 ) 


6.2 


0.0 


2.913 


230 


0.173 


0.124 


6.0 


0.0 


2.119 


219 


0.191 


0.133 


5.93 


0.0 


1.860 


214 


0.198 


0.136 


5.2 


0.1350 


1.876 


181 


0.213 


0.149 



Table 2: Coupling constants using A^g from QCDSF-UKQCD collaboration ]72^ |. The scale is set 
from tq = 0.5. 



in the light quark mass we don't use the mass of the D meson, as this is known to be 
problematic ]74l]. As we only have one sea quark mass value, the chiral extrapolation is 
only done on the masses of the valence quarks. Hence, we obtain a result for the partially 
quenched charm mass. 

In our earlier work [Q| we used the spin average of the pseudo-scalar and vector masses 
of heavy-light mesons, because the pseudo-scalar-vector mass splitting is underestimated 
in the quenched QCD. However, we found it difficult to reliably estimate the M2 meson 
mass for the vector meson, hence we now only use the pseudo-scalar meson. 

The meson masses are first interpolated to the strange quark mass using a simple linear 
ansatz 

M(mi) = ai + b l m l (5.1) 
where mi is the mass of the light quark. We denote this Mfj s - The strange quark mass has 



already been determined by the UKQCD collaboration for these ensembles [30, 75]. The 
masses are then interpolated to the charm mass using 

M(m Q ) = a h + b h m Q (5.2) 

where tjiq is the mass of the heavy quark. All the quark masses have corrections which scale 
with the quark mass, but the inverse hadron mass scales with quark mass, so it is unclear 
whether to plot arriQ vs. l/Mjjs or vs. Mh s q- We do the latter as we are interpolating 
in a finite range of am,Q, where a polynomial in arriQ can be expanded in terms of 1/amQ. 
Hence we are treating the charm quark as a heavy light quark rather than a light heavy 
quark (where HQET based extrapolations would be appropriate). We consistently include 
any mass dependent renormalization factors in the definition of tuq. Another option would 
have been to do the interpolation without the mass dependent renormalization factors, and 
then apply them after the fit. 

We use capital M for the meson masses and small m for the quark masses. In our 
analysis we consider the options: M\ versus mo, M\ versus m pcac , M\ versus mi, and M2 
versus m-2- The m-i mass has a different mass behaviour to that of m-2 in the heavy quark 
limit, so it only really makes sense to match the M2 mass with the m2 mass. 



6. Results 

In this section we include some examples of the data fits that were performed, present 
results at fixed lattice spacing and consider the systematic error associated with m c (m c ) 
in the continuum limit. 
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Figure 3: Example of a typical effective mass plot for the pseudo-scalar meson mass. The fit uses 
three exponentials to model data generated at (3 = 5.2, kjj = 0.1130, kl = 0.1340 and p = 1,0,0. 
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Figure 4: Example of a typical correlator plot for the PCAC mass. The fit uses a constant to 
model data generated at (3 — 5.2, kh = 0.1130, kl = 0.1340 and p = 0. 



6.1 Examples of data fits 

An effective mass plot for the pseudo-scalar and a plot of the PCAC correlators are pre- 
sented in figures H| and (^respectively. The heavy quark interpolation for vtla on the matched 
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ensembles is shown in figure || 



6.2 Results at fixed lattice spacing 

As noted in many places (for example |3(J), the determination of the M2 meson from the 
dispersion relation in equation 3.15 produces masses with bigger errors than for the M\ 
mass. So we use equation |3.2C to generate a perturbative estimate of the M2 mass. 

We test the approach in figures || and [7], where the perturbative estimate of M2 is 
compared against the non-perturbative estimate from the dispersion relation. At the finer 
lattice spacing, (3 = 6.2, all four definitions of the quark mass essentially agree. 



p 




m 


mi 


m 2 


5.2 
5.93 
6.0 
6.2 


1.327(4)i*j 
1.473(4)1^ 

1.439(4)1^ 

1.352(5)1|° 


0.952(1)1$ 
0.978(1)1| 
1.025(2)1^ 
1.147(3)±I 2 


1.247(3)lf 
1.253(3)1^ 
1.265(4)lf 
1.267(3)lf 


1.266(3)1^ 
1.274(2)1^ 
1.283(3)1| 
1.279(3)1* 6 



Table 3: The mass of the charm quark in the MS scheme at the charm mass scale, for different 
analysis techniques. We use ro = 0.5 fm as the central value and match at fi = 2/a. The first error 
is statistical and the second is due to the perturbative matching. 

Our results for the mass of the charm quark are presented in table ||[ The most striking 
point about the data is the differences between the results for the PCAC and vector masses. 
In quenched QCD the difference between the PCAC and vector quark masses is known to 
decrease as the continuum limit is taken |36|, 

We have also used the non-perturbative value for ca which was recently published 
in [42] with nf = 2 sea quarks. This reduced the lattice quark masses by 10%, consistent 



with the expectations in [42|. 

In figure || the dimensionless quantity r$M c (where M c is the RG invariant mass) is 
plotted against the square of the lattice spacing. We compare our data to that of other 
groups at non-zero lattice spacing. 

We used the convention for the renormalisation group invariant mass (torgi) used by 
the ALPHA collaboration |3£ 



ttirgi = rn(2b g 2 ' 



-d /2b 



exp 



1 



exp 



ch] 



1 1 

+ 



P(v) bo 7 ! 3 b\r] 



(6.1) 



This was different to the convention used in the RunDec package [FT]. 

The large splitting between the vector and PCAC definition of the quark mass in our 
unquenched data is seen to be consistent with the quenched data of Rolf and Sint. However, 
the use of the nonperturbative renormalization factors by Rolf and Sint || complicates the 
comparison. Figure || shows that the agreement between the final answer for the charm 
quark mass from Q and Rolf and Sint || is fortuitous. The final result from || was 
the average of the quark mass from the PCAC and vector currents. This is only a good 
estimate if the leading lattice spacing dependence from the PCAC and vector masses have 
opposite sign. 
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Figure 5: The hadron mass versus the bare quark mass for j3 = 5.2. The dotted line highlights 
the physical meson mass in lattice units. 



6.3 Estimating the final systematic error 

To set the lattice spacing we use ro between 0.5 to 0.55 fm. The advantage of using ro to 
determine the lattice spacing is that it is relatively easy to determine the value of ro in 
lattice units. In |^2[, the results for ro in physical units, from a number of different calcu- 
lations with unquenched Wilson like quarks were reviewed. All the results were between 
0.5 and 0.55 fm for ro- The new results from unquenched calculations using improved 



staggered fermions [76, 77, 78] are finding ro values of 0.462(H)(4) fm. The use of the 
HPQCD |78[ value for ro in the current generation of unquenched calculations with Wilson 
like fermions has been discussed in p3fl . The QCDSF collaboration have used the nucleon 
mass, in unquenched clover calculations, to estimate ro = 0.467 fm. However, it is diffi- 
cult to do a reliable chiral extrapolation of the nucleon mass in the current generation of 
dynamical lattice QCD calculations, that use clover fermions [jT9| 1, because of the size of 
the sea quark masses. For the data in this calculation we feel it is reasonable to use the 
estimate of ro between 0.5 and 0.55 fm. 

We now consider the continuum limit of the charm quark mass using the data in table |3|. 
One issue concerns the dependence of the mass of the charm quark upon the lattice spacing. 
One simple model for the dependence of the charm mass at non-zero lattice spacing is 



m c (m c ,a) = m c (m c )i + s/a 



(6.2) 
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Figure 6: The Mi and the M% mass computed from the dispersion relation, treel level and one 
loop in perturbation theory at (3 = 5.2. The horizontal line represents the physical mass in lattice 
units. 



The term linear in the lattice spacing comes from the use of improvement to one loop 
accuracy. We also tried a continuum extrapolation that was quadratic in the lattice spacing. 

m c (m c , a) = m c (m c ) q + s q a 2 (6.3) 

In table |I] we report the continuum extrapolation of the data in table |3[ In figures [| 



and 10 the charm mass is plotted as a function of lattice spacing with the fitted continuum 
extrapolation using the model in equation 3.S and |6.2| respectively. 



method 


mffrn,)! GeV 


si GeV 


m c (m c ) q GeV 


s q GeV 2 


Specie 


1.14(18) 


0.64(40) 


1.27(10) 


0.73(46) 


m 


1.45(5) 


-0.90(13) 


1.27(3) 


-1.0(2) 


mi 


1.29(7) 


-0.07(15) 


1.28(4) 


-0.08(17) 


m 2 


1.30(4) 


-0.05(8) 


1.29(2) 


-0.05(9) 



Table 4: Continuum limit of the mass of the charm quark mass in MS at the mass of charm for 
different analysis techniques. 

Figure || shows that a consistent continuum limit is obtained for all four definitions 



of quark mass if the extrapolations are done with equation 6.3, however it is difficult to 



give a rigorous argument in favour of this type of extrapolation. Figure 1C shows that the 
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Figure 7: The M\ and the Mi mass computed from the dispersion relation, tree level and one 
loop in perturbation theory at [3 = 6.2 quenched. The horizontal line represents the physical mass 
in lattice units. 



continuum limit of the PC AC and vector masses is inconsistent with the FNAL result, if 
equation |6.2| is used to take the continuum limit. This fit looks poor and we speculate that 
the continuum extrapolation should be done with a combination of linear and quadratic 
dependence on the lattice spacing. 

We also tried enforcing the same continuum limit for the vector and PCAC masses 
with fit parameters for both linear and quadratic terms in the lattice spacing. This gave 
VIS (m c ) = 1.57 ± 0.57 GeV. The fit is plotted in figure 10. Small O(a) terms are not 



rn 



obtained, as the fit finds that both the O(a) and 0(a 2 ) terms are large with opposite sign. 

The situation may have become clearer, if we had used quenched QCD calculations at 
finer lattice spacings. However this would just repeat the work of Rolf and Sint ||. To 
quote an unquenched result we need a formalism that produces the charm mass with a 
very weak dependence on the lattice spacing. As a check on our calculation, in table [5] we 
compare our results obtained with ALPHA formalism used by Rolf and Sint Q, with the 
results of Rolf and Sint. There is reasonable agreement between the two results. 

The continuum extrapolation of the m\ and 7712 masses is essentially consistent whether 
the numbers are extrapolated to the continuum limit either quadratically or linearly with 
lattice spacing. The extrapolation of heavy-light decay constants, obtained from calcula- 
tions that use the FNAL formalism, were also insensitive to whether a linear or quadratic 
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Figure 8: Recent data for the charm mass in the RGI scheme as a function of the lattice spacing 
The filled circles are the results for n/ — 2, all the other data is quenched. 



p 


This work 
m$ GI GeV m%j£ GeV 


Rolf and Sint 
m$ GI GeV m^gi GeV 


6.0 
6.2 


3.273(39) 4.430(52) 
3.768(45) 4.299(49) 


3.224(41) 4.331(59) 
3.711(47) 4.277(55) 



Table 5: A comparison of the results from this paper with those from Rolf and Sint j|. The 
renormalisation group invariant quark mass at charm for the vector and PCAC quark masses are 
reported using nonperturbative renormalisation factors. 



extrapolation in lattice spacing was done [|3C|]. 

For our quenched number we use the m\ number in the continuum limit from a linear 
extrapolation in lattice spacing. We also quote an error of 10% to account for variations 
in lattice spacing determinations. Hence, our final result for mass of the charm quark in 
the continuum limit of quenched QCD is 

m c (m c ) = 1.29(7)(13) GeV (6.4) 

The first error is a combination of statistics and an estimate of the error due to only using 
perturbation theory to one loop order. 

The data in table |3| do not show any clear pattern for the mass of the charm quark to 
differ between quenched and unquenched QCD at the fixed lattice spacing of 0.1 fm, with 
sea quarks with masses close to the strange value. Based on the scaling of the masses in 
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Figure 9: The quenched continuum limit of the scale invariant charm quark mass. The errors are 
from perturbative matching and statistics combined in quadrature. The curves are fits to a slope 
with 0(a 2 ) dependence only. 

quenched QCD, for the unquenched data we use the value of the mi mass with a 10% error 
for determining the lattice spacing. Our best nj = 2 number is 

m c (m c ) = 1.247(3)lf (120) GeV (6.5) 



7. Conclusions 

Our final result is m^ s (m c ) = 1.29(7)(13)GeV in quenched QCD. We found that the 
ALPHA and FNAL formulations gave consistent numbers, after making some assumptions 
about the lattice spacing dependence. Our result is consistent with previous results from 
quenched QCD given in table ||. 

We have determined the mass of the charm quark in two flavour QCD at a lattice 
spacing of 0.1 fm with a sea quark mass around strange value. We did not observe any 
unquenching errors. 

The most important task for future unquenched lattice calculations is to control the 
lattice spacing errors in the mass of the charm quark. The large scaling violations found in 
the charm quark mass (see figure |8|) with the improved clover action in quenched QCD can 
be controlled by the brute force method of using a lattice spacing of 0.05 fm. The timing 
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Figure 10: The quenched continuum limit of the scale invariant charm quark mass. The errors are 
from perturbative matching and statistics combined in quadrature. The curves are fits to a slope 
with 0(a) dependence only. The curve labelled by pcac + mo is a fit where the pcac and vector 
quark masses are forced to have the same continuum limit. 



Group 



m 



MS 



(m c )GeV 



estimates in [64, 35, 81] suggest that this approach will be not be easy for unquenched 
calculations because of the large computational cost in reducing the lattice spacing. 

The FNAL formalism seems to have a 
better scaling behaviour than the ALPHA 
formalism for this data set. We note that one 
of the PCAC quark definitions used by Rolf 
and Sint also has a very weak lattice spacing 
dependence ||. The development of better 
fermion actions for heavy quark calculations 
is clearly desirable 



Another important systematic error that 
must be reduced originates from matching 
the lattice renormalisation scheme onto the 
continuum 



This work 
Becirevic et al. 
Rolf and Sint j 
Juge @ 
Kronfeld §§ 
Hornbostel et al. 
de Divitiis et al. 



1.28(3)(13) 
1.26(4)(12) 
1.301(34) 
1.27(5) 
1.33(8) 
1.20(4)(11)(2) 
1.319(28) 



Table 6: Mass of the charm quark mass from 
various quenched lattice QCD calculations. 



As we discussed in section |4,2j , 
there are many lattice techniques for reducing the error on the matching of the quark 
masses to the MS scheme. Many of these techniques will benefit from unquenched data 
with finer lattice spacings. 
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